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Abstract 

The depth profile of the intrinsic magnetic properties in an Fe/Sm-Co bilayer fabricated under 
nearly optimal spring-magnet conditions was determined by complementary studies of polarized 
neutron reflectometry and micromagnetic simulations. We found that at the Fe/Sm-Co interface 
the magnetic properties change gradually at the length scale of 8 nm. In this intermixed interfacial 
region, the saturation magnetization and magnetic anisotropy are lower and the exchange stiffness 
is higher than values estimated from the model based on a mixture of Fe and Sm-Co phases. 
Therefore, the intermixed interface yields superior exchange coupling between the Fe and Sm-Co 
layers, but at the cost of average magnetization. 

PACS numbers: 75.70.Cn, 61.05.fj 
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I. INTRODUCTION 



Exchange-coupled, high-magnetization (soft) and high-anisotropy (hard) magnetic phases 
have potential applications as both ultra-strong permanent magnelP and ultra-high-density 
recording mediaP While the intrinsic properties of the two phases play the most impor- 
tant role, optimization of the interface properties are also important to achieve the best 
performance. For example, the interface morphology need to be optimized for good ex- 
change coupling between the soft and hard phases. Micromagnetic simulations suggested 
that magnetically graded interfaces, whose magnetic properties are gradually changed over a 
distance ~ 10 nm, is more effective than sharp interfaces, in order to increase the nucleation 
field in the soft phase and decrease the switching barrier of the hard phase.l^ Magneti- 
cally graded interfaces can be fabricated from chemically intermixing phases via nanotech- 
nology.l^ Nanoscale spatial resolution is needed to experimentally determine the intrinsic 
magnetic properties of hard-soft heterostructures. Interestingly, for Fe/Sm-Co spring mag- 
nets under optimized fabrication conditions for the maximum energy product, the Fe layer 
and the Sm-Co layer have an intermixed interface over a length scale ~ 5 — 10 nm.^ ^ Al- 
though the intrinsic magnetic properties of individual Fe films and Sm-Co films have been 
well studied,'^ there lacks quantitative knowledge about the intrinsic magnetic properties of 
the intermixed Fe/Sm-Co interface due to its complex composition.'^ 

Polarized neutron reflectometry (PNR) is a sensitive tool to study the depth profile of 
magnetic structures within multilayers with sub-nanometer depth resolution,'^^'^ and has 
been used to study exchange coupling and anisotropy.'^^'^ For example, in a NiFe/FePt 
spring magnet, O 'Donovan et al. have determined that the twisted magnetic structure of 
the spring magnet is not confined to the magnetically soft layer, but also penetrates into 
the hard magnetic phase.'^ Recently, Kirby et al. showed qualitatively that the structural 
gradation yields a graded anisotropy in Co/Pd multilayers.'^ 

In the present work, we report PNR studies on a Fe/Sm-Co bilayer, which was fabricated 
under nearly optimized spring magnet conditions. We first confirmed that there is a struc- 
turally intermixed interface of ~ 8-nm wide between the Fe and Sm-Co layers via combined 
X-ray reflectometry (XRR) and PNR studies. We also determined the depth profile of the 
saturation magnetization Ms- Furthermore, the profiles of the magnetic anisotropy K and 
the exchange stiffness A were obtained by analyzing the PNR data with the aid of micro- 
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magnetic simulations. The magnetic properties in the intermixed region were compared to a 
model based on a mixture of Fe and Sm-Co phases. The saturation magnetization is slightly 
lower than the value estimated from the model, suggesting new compounds formed in the 
intermixed region. The intrinsic anisotropy is also lower than the value from the model. 
However, the exchange stiffness is higher, so that the interface efficiently couples the Fe and 
Sm-Co layers. 

II. EXPERIMENTAL METHODS 

The Fe/Sm-Co thin film was fabricated via dc magnetron sputtering onto an MgO (110) 
single-crystal substrate with a nominal structure of (10 nm Cr)/(10 nm Fe)/(20 nm Sm- 
Co)/(20 nm Cr)/MgO. The sample size is 10 x 10 mm^. The Cr (211) buffer layer and the 
Sm-Co layer were grown at 400 °C, and the Fe layer was grown at 100 °C.'^ The Sm-Co 
layer, nominally Sm2Co7 composition, has an uniaxial, in-plane magnetic easy axis.l^ Mag- 
netic hysteresis loops were obtained by means of vibrating sample magnetometer (VSM). 
XRR studies were performed with a X-ray diffractometer using Cu radiation. PNR 
experiments were conducted at the SNS at Oak Ridge National Laboratory, at the Mag- 
netism Reflectometer.'^ This is a time-of-flight (TOF) instrument with a wavelength band 
2 - 5 A and the polarization efficiency of ~ 98%. All experiments were conducted at room 
temperature. 
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FIG. 1. (Color online) Easy-axis magnetic hysteresis loop. The film is satmated at 1 T. The two 
circles label the magnetization states measured by PNR, with applied fields of +1.11 and +0.39 T, 
respectively. 

The external field was applied along the magnetic easy-axis for PNR experiments. For 
convenience, we define the field direction as +x and the sample surface normal as —z. 
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Figure [T] shows the easy- axis magnetic hysteresis loop. The film saturates at fields above 
+1 T. There are two major reversal processes along the increasing field branch, centered 
at +0.35 and +0.75 T, respectively. A minor soft-phase was also observed, which reverses 
below +0.05 T. PNR measurements were performed at +1.11 T and in a demagnetized state 
at +0.39 T, after saturation with a —1.11 T field, in order to determine the depth profiles 
of the saturation magnetization, the exchange stiffness and the magnetic anisotropy. 

Reflectometry is a non-destructive method to determine scattering length density (SLD) 
profiles. In the specular condition, XRR yields the depth profile of the electron density, 
which can be used to reconstruct the chemical structures. Neutrons interact with both nu- 
clei and the internal magnetic field. Spin-up and spin-down neutrons feel the same nuclear 
scattering potential, but an opposite magnetic scattering potential. From subsequent mea- 
surement with oppositely polarized neutron beams, the two contributions can be separated 
in order to reconstruct both depth profiles of the chemical structure and the magnetization 
vector .fim^lxhe contrast of the neutron SLDs between Fe and Sm-Co is high, so that the PNR 
experiments are sensitive to the interfacial structure of interest. There are four refiectivities 
in PNR, including two non-spin-fiip (NSF) refiectivities R'^'^ and R , and two spin-fiip 
(SF) refiectivities R^~ and R~~^. SF scattering occurs when the sample's magnetization 
vector has a non-zero component perpendicular to both the neutron's polarization direction 
and the momentum transfer direction. During magnetization reversal, the SF scattering 
occurred in a sufficiently high magnetic field so that it was necessary to take into account 
the Zeeman effect when analyzing the data.'^ (See Appendix [A| ) 

For PNR, Fredrikze's formalisnP^ was used to determine the optics' polarization effi- 
ciency and the direct beam (DB) spectra from the measured four reflectivities of the DB. 
The polarization corrections with the error propagation were made following Wilder's for- 
malism.'I^ Simulations of the reflectivities were based on the Parratt formalism.®! A rough 
interface was modeled as a sequence of very thin slices whose SLDs vary, following an error 
function so as to interpolate between adjacent layers.'^ The instrumental resolution was 
handled by Gaussian convolution. The GenCurveflt programpH using the genetic algorithm 
was employed for the model optimization. 
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III. EXPERIMENTAL RESULTS AND ANALYSIS 
A. Depth profile of Saturated Magnetization 




FIG. 2. (Color online) (a) XRR data and (b) PNR data taken at saturation field +1.11 T. 
(c) Depth profiles of neutron nuclear SLD (thin black line, left), neutron magnetic SLD (thick red 
line, left) and X-ray SLD (dashed green line, right). Only the real parts of the X-ray and the 
nuclear SLDs are shown. 

The XRR data and the specular PNR data taken at a saturation field of +1.11 T are 
displayed in Figs. [2^ and |2]d, respectively. The reflectivities are shown as functions of the 
momentum transfer perpendicular to the film plane Q-, = AnsinO/X, where 6 is the incident 
angle and A is the wavelength of the radiation source. Since there is no spin flip scattering 
exists at saturation, we measured only {=R~^~^) and R" {=R ). Because of the sample's 
high saturation magnetization, R'^ ^ R~ for Qz above the critical edges {Qz > 0.018 A~^). 
There is an interesting feature in the PNR data below the critical edges: both R'^ and 
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R~ showed frustrated total reflection. This is due to the enhanced absorption at resonant 
conditions.'^ Simulations show that the Q separation between the two dips strongly depends 
on the magnetization along the field direction, especially that of the Sm-Co layer, which is 
the locus of absorption. Therefore, this separation, as well as the splitting between R'^ and 
R~ above the critical edge, are direct indications of the sample's magnetization. 

The chemical and the magnetic SLD profiles were determined by fitting the XRR and the 
PNR data simultaneously with the same structural model.l^ Since there are well tabulated 
values for the neutron and the X-ray scattering length and absorption length for each ele- 
ment, both the real and imaginary parts of the X-ray and the neutron's nuclear SLDs were 
calculated from the chemical composition and the film density. Therefore a single parameter 
was used to determine both the X-ray and the neutron nuclear SLDs for each layer that has 
a well defined composition. This parameter is the mass density for that particular layer. 
However, the intermixed interface has a complex chemical composition that extends > 5 nm 
in depth.'^ The X-ray and the neutron's nuclear SLDs in the intermixed region need to be 
determined separately, which was done by introducing an intermixed layer between the Fe 
and Sm-Co layers. Therefore, both the nuclear and the magnetic SLD profiles can vary 
in more sophisticated ways and they are less correlated to each other over the intermixed 
region. Actually, the peak and dip positions in R~ between 0.025 and 0.045 A can not 
be reproduced without introducing this intermixed layer, which indicates that a single er- 
ror function is not sufficient to model the nuclear and/or the magnetic SLD profiles at the 
Fe/Sm-Co interface. 

The best-fit curves for the reflectivity data are overlaid on the data, and the SLD depth 
profiles are shown in Fig. The SLD profiles show sharp transitions at the Cr/MgO and 
the Sm-Co/Cr interfaces with rms roughnesses of ~ 0.4 nm for both interfaces. However, the 
SLD profiles gradually change between the Fe and Sm-Co layers over a distance of ~ 8 nm. 
Such a large intermixed region is consistent with previous energy-dispersive-spectroscopy 
(EDS) results from samples fabricated under the same conditions.'^ In the intermixed region, 
there appears a shoulder in the nuclear SLD profile (Fig. showing that it is chemically 
rich in Fe. However, this feature is not present in the magnetic SLD profile. This will be 
discussed below. The parameters of interest are listed in Table |T} which shows that both the 
Fe and the Sm-Co layers have their saturation magnetizations close to the literature valued 
despite the large intermixing at the interface. 
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TABLE I. The layer thickness, rms interface roughness, the nuclear and magnetic SLDs, saturated 
magnetization Ms, exchange stiffness A and uniaxial anisotropy K in the layers of interest from 



the best fits. Nuclear and magnetic SLDs and Ms are from fitting the +1.11-T data (Sec. Ill A). 



A and K are from fitting the +0.39-T data (Sec. IIIB). Typical literature values are Ms = 
1700 (550) emu/cc, A = 2.8 (1.2) ergs/cm, and K = 10^ (5 x 10^) ergs/cc for the Fe (Sm-Co) 
layer.l^ (*Literature value of Kpe was used during the optimization.) 





Fe 


mixed 


Sm-Co 


Thickness (nm) 


7.4 ±0.3 


4.2 ±0.3 


20.1 ±0.2 


Roughness (nm) 




4.3 ±0.3 2.6 ±0.2 




Pnuc (10-6 A-2) 


7.9 ±0.2 


5.2 ±0.2 


1.64 ±0.04 


Pmag (10"^ A-^) 


4.9 ±0.1 


2.5 ±0.1 


1.62 ±0.03 


Ms (emu/cc) 


1700 ± 50 


890 ± 50 


570 ± 10 


M+0.39 T (emu/cc) 


1540 ± 10 


730 ± 30 


420 ± 10 


A (10~6 ergs/cm) 


2.6 ±0.1 


2.6 ±0.1 


1.2 ±0.1 


K {10^ ergs/cc) 


0.001* 


3.0 ±0.2 


27 ± 1 



B. Depth Profiles of Exchange Stiffness and Uniaxial Anisotropy 

In order to get insight into the intrinsic magnetic properties, besides Ms, in the intermixed 
region, we also studied the magnetization structure in a demagnetized state. The PNR data 
were collected at ±0.39 T after the negative saturation in —1.11 T. The field is sufficiently 
high so that the specular SF reflection showed at off-specular positions on the detector due 
to the Zeeman effect, as shown in Fig. [6] Therefore, in the analysis the momentum transfer 
Qz was modified according to Eq. |Alf ^^ * ^^ l It is worth noting that the off-specular scattering 
in the SF reflectivity is dominated by the Zeeman effect in our data and possible off-specular 
scattering due to in-plane magnetic domains is not distinguishable from the background, and 
is therefore not considered in the data. 

The PNR data are shown in Fig. [3^. In contrast to the saturation case, R is higher 
than W"^ for most Qz values above the critical edge, which indicates that the average 
is still along the negative applied field direction. Both the splitting between R and 

7 



R^^ and the separation between the dips below the critical edge are smaller because the 
average is much lower than Mg. At the same time, there is significant SF scattering 
of the same amplitude as the NSF scattering, indicating a large in-plane magnetization 
component perpendicular to the field direction (My). This is caused by the magnetization 
spiral structure observed in exchange-coupled bilayers during demagnetization.'^ Figure 
plots the depth profile of the magnetization vector M that yields the best fit, which is 
obtained from micromagnetic simulations. Our PNR studies do not reveal the chirality of 
M so that a positive My is used for convenience. As expected, the magnetic moment of the 
Sm-Co layer is along the negative field direction and the magnetization vector rotates into 
the intermixed region and the Fe layer. 

Rohlsberger et al. have shown that the depth profile of the magnetization vector in 




FIG. 3. (Color online) (a) PNR data taken at +0.39 T after saturation in —1.11 T. W and 
R ^ are offset by a factor of 0.1 for clarity, (b) The depth profiles of the magnetization vector 
M (obtained from micromagnetic simulations) that yields the best fit to the PNR data. The field 
direction is along +x. 
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exchange coupled bilayers agrees with the micromagnetic simulation from the ID spin-chain 
model.'^Fitzsimmons et al. further determined the micromagnetic parameters for individual 
layers in the exchange coupled DyFe22/YFe2 superlattice using combined studies of PNR 
and micromagnetic simulations.^^ For the Fe/Sm-Co spring magnet studied here, the Fe 
and the Sm-Co layers are largely intermixed for ~ 8 nm at the interface and it is naturally 
to expect that the intrinsic magnetic properties, i.e. uniaxial anisotropy K and exchange 
stiffness A, vary gradually due to the gradation of the chemical composition. Hence, we 
extended Fitzsimmons' approach and constructed depth profiles of the intrinsic magnetic 
properties. Given a profile of the micromagnetic parameters, the equilibrium magnetic 
structure can be simulated by energy minimization via the ID spin-chain model,'^ from 
which the neutron's magnetic scattering potential can be directly calculated followed by the 
spin-dependent neutron reflectivities. Therefore, the intrinsic magnetic properties can be 
optimized to yield the best fit to the experimental PNR data. 

These profiles of K and A are built up of layers identical to those of the chemical structure 
as determined from the fits to the XRR and the PNR data in saturations. The gradual 
changes between the layers are computed using the interface roughness in the same way as 
is done for the SLDs. When fitting, the chemical structure and the roughnesses are fixed to 
those already determined. M is allowed to be lower than Mg to account for the effect from 
the minor soft phase mentioned above. The model optimization is not sensitive to the i^^e 




10 20 30 40 

z (nm) 

FIG. 4. (Color online) Top: The depth profile of the nuclear SLD in gray scale in the unit of 
10^^ A^^. Bottom: The depth profile of the exchange stiffness A (left, red solid line) and the 
uniaxial anisotropy K (right, blue dashed line) which gives the best fit to the PNR data. 
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because the associated magnetic energy is negligible, and therefore the literature value of Kpe 
is used. The error bars of the parameters only reflect the statistical error and the accuracy 
of the values depends on the model. Therefore it is important here to have the additional 
layer between the Fe and the Sm-Co layers in the model because the exact mapping function 
from the chemical structure to the micromagnetic properties is unknown. The layer has its 
own free micromagnetic parameters so that M, A and K in the intermixed region are able 
to vary more independently from those in the Fe layer and the Sm-Co layer. We checked the 
robustness of the model by allowing different interface roughnesses for M, A, and K and 
found that the optimized micromagnetic parameters have similar depth dependencies. 

The depth profiles of A and K from the best fit are shown in Fig. |4] and the values of 
M, A and K for layers of interest are listed in Table |T} Both the parameters A and K show 
monotonic depth dependence as expected. Msm-Co is ~ 20% lower than the saturation 
value, from which we estimated that ~ 10% of the Sm-Co was reversed due to the minor 
phase. Therefore, the ID spin-chain model is approximately valid and yields fairly good fits, 
which were overlaid on the data in Fig. [3^. As also listed in Tabled in the intermixed region, 
A is close to that of the Fe layer, but K is much lower than the average value between the 
Fe and the Sm-Co layers. 

IV. DISCUSSION AND SUMMARY 

By assuming that the intermixed region is composed of a mixture of the Fe and Sm-Co 
phases, the saturation magnetization can be computed from the relative volume of both 
phases via M = fsMg + fhMh, where fs [Ms) and fh (Mh) are the relative volume (the 
saturation magnetization) of the Fe phase and the Sm-Co phase, respectively. From the 
depth profile of the nuclear SLD, the relative volume of the two phases can be determined. 
Following this approach, the calculated saturation magnetization profile is shown in Fig. [5^. 
Also shown is the saturation magnetization profile calculated from the magnetic SLD profile. 
Clearly, the mixture model predicted slightly larger saturation magnetization in the inter- 
mixed region. This agrees with the suggestion that new compound(s) is (are) formed in the 
intermixed region,'^' possibly due to interdiffusion between the Fe and Sm-Co layers. The es- 
timation shows that the average magnetization of the whole sample is reduced by ~ 4% due 
to the intermixing of the two components. The exchange stiffness and the anisotropy of a two 
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FIG. 5. (Color online) Depth profiles of (a) the saturation magnetization, and (b) the exchange 
stiffness (left) and the uniaxial anisotropy (right) calculated by assuming that the intermixed region 
is a mixture of the Fe phase and the Sm-Co phase (dashed lines). The experimentally determined 
profiles are plotted as solid lines. Also shown in the top panel is the depth profile of the nuclear 
SLD in gray scale. The unit is 10~^ A~^. 

phase mixture can also be approximately estimated by the relative volume: K = fgKg + fhKh 
and A = f,A, + UAhW The results are shown in Fig. [5|d. In comparison to those deter- 
mined from the PNR experiments, the model predicts higher anisotropy but lower exchange 
stiffness in the intermixed region. Therefore, we conclude that the optimized Fe/Sm-Co 
interface for a spring magnet reduces the average magnetization. 

Interfaces with graded magnetic properties are desirable for superior performances in ex- 
change coupling systems, such as spring magnets and magnetic media.'^ A straightforward 
approach to achieve magnetically graded interfaces is to make structurally graded inter- 
faces.lSHZI The micromagnetic properties could be predicted from the composition. However, 
it is challenging if multiple phases coexist even when the micromagnetic properties of each 
individual phase are known, due to lacking of prior knowledge of the interphase coupling. 
As shown above, the intrinsic magnetic properties in the intermixed Fe/Sm-Co interface 
deviates from the results predicted by the mixture model. Therefore, a complex correlation 
between the magnetic properties and the nominal composition may exist at the chemically 
graded interfaces. 
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In summary, we determined the intrinsic magnetic properties of a Fe/Sm-Co bilayer un- 
der the nearly optimized fabrication condition for spring magnets. We determined that the 
intermixed interface between the Fe and Sm-Co layers extends ~ 8 nm, where the magnetic 
properties changes gradually, as expected. We compared the magnetic properties with the 
prediction of a mixture model and found, in the intermixed region, the saturation magne- 
tization is slightly lower than that estimated from the model, but the exchange stiffness is 
higher. This observation indicates that the intermixed interface is efficient for magnetically 
coupling the Fe and Sm-Co layers but at the cost of the average magnetization. The intrin- 
sic anisotropy is also lower than the value from the model. Overall, the intrinsic magnetic 
properties in the structurally intermixed region may not be predicted correctly by the mix- 
ture model of nominal compositions, which is worth keeping in mind when designing the 
magnetically graded interfaces. 
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Appendix A: Zeeman effect and ofF-specular scattering 

The Zeeman effect in PNR was firstly clarified by Felcher et. alJ^^ The effect shows 
up in the SF refiectivities when the difference of the Zeeman energy for spin-up and spin- 
down neutrons is not negligible in a sufficiently high magnetic field. Let k± and k// label 
the perpendicular and the parallel components of the wavevector with respect to the sample 
surface in the vacuum region, i.e., before neutrons enter the magnetic-field region of interest. 
k± = 2'Ksin9/\ ^ k/j at glancing angles. There is negligible energy exchange between the 
sample and the neutrons for elastic scattering, therefore the Zeeman energy change after SF 
scattering accompanies a change of the kinetic energy, which is associated with k. The energy 
change after spin-fiip is twice the Zeeman energy. If the sample is magnetically homogenous 
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in the film plane as seen by the neutrons, the energy change is totally associated with k±. 
Since k±s are small quantities, the outgoing angles of the spin- flipped neutrons are different 
from the incoming angles. Therefore the SF reflections appear at off-specular positions, 
which have a characteristic field dependence .'^^'^ The Zeeman effect on the SF reflections 
becomes clear in a magnetic field on ~ 0.1 T or higher. Figure |6] shows an example. It 
is clear that the locus of the SF reflections deviates away from the specular position, but 
are described by the prediction after considering the Zeeman effect. Beside the reflection 
due to the Zeeman effect, there is no noticeable off-specular scattering in the SF channels; 
therefore, the scattering from the in-plane magnetic domains are not considered during the 
data analysis. Consequently, the momentum transfer Qz are no longer equal to 2k± for 
spin-flipped neutrons, but follows 

Qt- = + Vkl + 2CB ^^^^ 
Q-+ = k± + ^kl - 2CB, 

where C = |2m„/in//i^| = 2.906 x 10^^ A^^/T and B = hqH is the applied magnetic field. 
R ^ is forbidden when k± < \/2CB because there is no enough kinetic energy to compensate 
the Zeeman energy change.!^ At the same time, the minimum Q^" is also \/2CB for finite 
reflectivity. Therefore, there is a cutoff for non-zero SF scattering, QsFcutoff = V2CB. 

Actually, the kinetic energy also changes when neutrons enter and leave the magnetic-field 
region of interest, but the change of k is essentially in k// in these cases, which changes 
the time-of-flight. However, the relative change of k// is so small in laboratory fields, and 
TOF only changes ~ 10 ns, while the width of the neutron pulse is ~ 0.2 ms. Therefore, it 
does not result in any observable affect. 

SF reflection is typically weak when the Zeeman effect matters, therefore the Zeeman 
effect is usually not considered during data reduction and analy sis .H'^ However, we observed 
significant SF scattering at +0.39 T, which is sufficiently high so that the Zeeman effect needs 
to be considered. We adopted the generalized algorithm of the GEPOREJ^ for the PNR 
simulations. B rather than M is used to calculate the magnetic scattering potential since 
H is now comparable to M and the spin eigenstates are aligned along B rather than M. 
Considering the problem along the normal direction of the sample surface (let it be z), it is 
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FIG. 6. (Color online) The contour map of the nentron intensities measured with the polarization 
analysis, which are presented as functions of the neutron wavelengths A. The data were collected 
from a Fe/Sm-Co spring magnet sample at an incident angle 9i = 0.24° in an external field of 
+0.39 T after a negative saturation. The dashed lines indicate the specular reflection positions 
without considering the Zeeman effect, i.e. 6f = 9i. The solid lines in the SF channels indicate the 
off-specular scattering following the the Zeeman effect, i.e. sin^ ^/ = sin^ 6i ± ^p^. 



reduced to a pair of coupled ID differential wave equations 



-^^2 -L V 



[z)-E 




0. 



(A2) 



vI/_(z) + y_+M/+(z) 

V is the potential operator. '^±{zys are the wave functions for spin- up and spin down states. 
/cq is the wavevector in vacuum, which can be either a real number or a pure imaginary 
number. States with k^^ ± Ci? < correspond to evanescent waves outside the sample. 
Due to the Zeeman energy, the coupled eigenstates of and have different wavevectors 
outside the sample, namely a/^q 2 - CB and + CB, which has been considered to set 
up the boundary conditions at the sample surface in the modified algorithm to take account 
the Zeeman effect. 
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